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ABSTRACT . 


This report describes a computer program which calculates three- 
dimensional viscous effects on the lift and pressure distributions of 
arbitrary three-dimensional bodies. The program is a combination of a 
panel method, which computes the potential flow about arbitrary three- 
dimensional lifting configurations, and a three-dimensional boundary-layer 
method, which calculates the viscous effects with small crossflow. These 
effects are applied to the three-dimensional shape, as in Part 1, ina 
"strip-theory" sense and the resultant viscous lift and pressure distribu- 
tions are again produced. The method of simulating the boundary layer in 
the final potential-flow calculation is by the addition of the displacement 
thickness to the original shape. 


The computer program is written in Fortran IV for the IBM 370 systems. 
16 temporary external units are used for storage. The region size needed 
to execute the program is about 360K bytes, but this is a direct function 
of the number of elements defining the configuration. 


Also presented in this report is a detailed description of the program 
logic, complete instructions for executing the program, and a sample case. 
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1.0 INTRODUCTION 


This is the second part of the user's manual for operating the combined 
three-dimensional potential-flow boundary-layer computer program. This part 
of the report concentrates on the procedure for obtaining the displacement 
thickness from the boundary-layer program in a three-dimensional form which 
allows for small crossflows. 


The operation of this program is similar to that described in Part 1 [2]. 
The flow field data are generated and the input body shape is separated into 
lower and upper regions. The flow field data are then input to the stream- 
line program to obtain the streamline curvature, orthogonal curvature and the 
streamline distances between the streamline points. These quantities are 
saved and reorganized to the form that is acceptable to the three-dimensional 
boundary-layer program which is based on that of reference 4. 


When the boundary-layer properties with small crossflows have been 
calculated, the displacement thicknesses, 6*, associated with the streamline 
points are obtained. The lower and upper surface 6* fields are defined at 
the streamline points. A surface interpolation method is used to obtain the 
6* at the control points that are used to compute the potential-flow calcula- 
tion. After the interpolated §* procedure is accomplished, the surface 
displacement method of boundary-layer simulation is used as described in 
Part 1 of this manual. The generated coordinates are the coordinates with 
the boundary-layer effect. These new body points are once again input to the 
potential-flow program to obtain the final potential-flow solution. 


At the present time, a multistep computer execution procedure is used 
for the above computation. This gives the user the flexibility to run the 
whole procedure as a single computer job, or to run each program step sepa- 
rately and examine the output at the end of each program step; in the latter 
case, data are saved and the user has the option of using it as input to the 
next step or to override it by card input. 


All the programs should be compiled and stored before executing with 
data. The programs are written in Fortran IV language to be run in the 
IBM/370 system and a direct access unit is required in the streamline 
calculation. 


2.0 PROGRAM NOTES 


2.1 Program Logic Flow Diagram 


(1) POTENTIAL FLOW CALCULATION 
(2) PREPARE FLOW FIELD DATA FOR 
STREAMLINE CA ATION 


STOP 


(1) STREAMLINE CALCULATION 
(2) PREPARE DATA FOR SMALL CROSS- 
FLOW B.L. CALCULATION 


STOP 


(1) SMALL CROSSFLOW B.L. 
COMPUTATION 

(2) GENERATE 6* FOR EACH STREAM- 

LINE POINTS 


STOP 


(1) REORGANIZE THE 6* DATA 
(2) BACK INTERPOLATION AT THE 
CENTROID LOCATION 


TERROR 


#0 STOP 


USING THE NEW 6* TO OBTAIN NEW 
COORDINATE 
SET _BL<0 


STEP 1 


STEP 2 


SEP: 3 


STEP 4 


2.2 Potential-Flow Calculation 


The potential-flow calculation is performed initially to generate the 
flow-field data, i.e., the surface velocities, pressure coefficients, centroids 
and unit normals, etc. This calculation is used again after the displace- 
ment thickness has been added to the input body coordinates. 


The major part of the potential-flow program is the same as described in 
Part 1 of this manual. There are a few additional routines here designed 
solely for the small crossflow version which are not needed in the version 
described in Part 1. 


The functions of these particular routines are: 


(1) Compute the u, v coordinate of each input element and their 
derivatives with respect to x and y. 


(2) Compute the stagnation point and separate the input lifting body 
into two regions. Each region starts at three elements before the 
stagnation point and goes toward the trailing edge. 


(3) Organize the flow-field data of these regions to be used in the 
streamline routines. 

The following units are saved: 

(1) Unit #2 — stores the u,v data. 


(2) Unit #21 — stores the flow-field data for each region. 


2.3 Region Data Organization 


A lifting section is separated into two reqions; namely, the lower surface 
region and the upper surface region. Each region starts ata location near 
the stagnation point and goes toward the trailing edge. A lifting strip, there- 
fore, is divided into two parts and a region is composed of the parts of these 
lifting strips. This is illustrated in the following sketch. 
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For each element in the region, there are 21 flow quantities associated 
with it, which over the whole region define the flow field that will be used 
for surface interpolation. The flow data and their storing sequence for 


each element follows. 


ee 


GIDE 
c 
(2) ve centroid component of the element. 
(3) 2, 
(4) MU 
(5) MV equations (9) and (15) of [1]. 
Ve | velocities parallel to & and n axes. See 9.3 of [1]. 
a Nn 
(8) Ms 
(9) vy components of surface velocity 
(10) ¥ 
He Uy | nonorthogonal coordinates pertaining to the midpoint of 
13 a panel. 
t) 
(13) an 
(14) ay 
(15) a13 
(16) any 
(17) ano components of the transformation matrix 
(18) an3 
(19) a3] 
(20) a5 
(21) a33 


These flow data are stored for each on-body element,strip by strip, from 
the starting element of the region to the trailing-edge element (excluding 
wakes). The lower surface is considered as regicn No. 1 and the upper surface 
is region No. 2. The program uses unit No. 21 to store the flow-field data. 
In addition to the flow-field data, eight quantities that are associated with 
the stagnation point for each strip are also stored following the flow-field 
data. They are: 


(1) CPSTAG pressure coefficient 

(2) WESTAG surface velocity 

(3) XSTAG 

(4) YSTAG components of the stagnation point 
(5) ZSTAG 
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(6) TSX 
(7) TSY components of the tangent vector 
(8) TSZ 


The aforementioned data are all to be used in the streamline routine. 
Specifically there are three external units to be saved: 


(1) Unit #4 — geometry data of each element 
(2) Unit #2 — nonorthogonal coordinates 
(3) Unit #21 — flow-field data of the two regions 


2.4 Streamline Calculation 


The streamline calculation program is a modified version of the streamline 
calculation from the Mark IV supersonic-hypersonic arbitrary-body program [3]. 
A direct-access unit is used to store the flow-field data (the unit is defined 
as unit #10 in the program). Pointers for regions, and streamline points are 
saved so that users may use the same data again if the streamline program is 
to be run independently. The flow field data can be input by cards or passed 
from the potential-flow calculation program in unit #21 as stated previously. 


After the flow-field data are input, an option is provided to check 
the flow-field data without going into the streamline calculation. If the 
streamline calculation is desired, the spacing between streamline points has 
to be input. The program will compute the streamline points from the trailing 
edge element toward the stagnation point using > fourth-order Runge-Kutta 
integration method with the user supplied spacing. Surface interpolation is 
performed to obtain the interpolated flow-field data at the streamline points. 
Orthogonal curvature is also computed at the same time. 


Upon completion of all the streamline points for both regions, the program 
will organize the streamline data such that the order is from the stagnation 


point to the trailing edge. Streamline curvatures are also computed. 


The following sketch illustrates the way the streamlines are calculated 
and organized for input to the boundary-layer calculation. 
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2.5 Boundary Laver with Small Crossflow Ca culation 

The small crossflow boundary-layer program used in this calculation is a 

modified version of [4]. The boundary-layer program uses the generated output 
from the streamline program, one streamline at a time, to compute the dis- 
placement thickness at each of the streamline points. 


The input to this proaram requires the following data: 


(1) Number of points, station (x), the arc length between these stations, 
their curvature and their velocities. 


(2) Reynolds number per foot, reference velocity in ft/sec. 


(3) The point number of the transition location on each strip. 


Items (1) and (2) above may be input by card as an option (see Section 2.3), 
as well as passed by the program. Only the transition location numbers must 


be input as punched cards. 


On output, the displacement thicknesses (6*) are stored for each stream- 


line point to be used for back interpolation. 


2.6 Back Interpolation 

The term back interpolation used here means that the displacement thickness 
values obtained from the small crossflow boundary-layer program at the streamline 
points, which define the 6* over the whole region surface, are used to 
obtain 6* at the original potential-flow calculation control points by using 


the surface interpolation method as illustrated in the following figure. 


STREAMLINES 


POTENTIAL FLOW 
CONTROL POINT 


After the <6* at each control point is obtained, they are added to 
the original input coordinates in the same manner as described in Part I, [2] 
and in [1]. 


2.7 Simulation Program 


This part of the program handles the generation of new coordinates with 
the 6* added on as described in Part 1 [2]. Output may be saved on external 
units or punched on cards. They are used as input when recalling the potential- 
flow program, 


2.8 Program Option 


Listed below are the options that may be selected. 


(1) All of the potential-flow calculation options (note the restrictions 
mentioned in the next section). 


(2) Options to make a complete run in the streamline program or partial 
run to examine the data. Option to print every streamline point 
or every other point. 


(3) Option to input the flow-field data by cards or by external unit, 
i.e., tape or disk. 


(4) Print output option in the small crossflow boundary-layer program. 
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CS, 


(5) 


(6) 


Option to execute the individual streamline in the small crossflow 
boundary-layer program. 


Option to override the internally generated transition location. 


Program Restriction 


At present, this program is subject to the following restrictions: 


(1) 
(2) 
(3) 
(4) 


(5) 


Applicable only to the lifting body. 
One angle of attack. 
One section input. 


The number of elements that define a region must not be greater 
than 500. 


No more than 100 points per streamline. 


3.0 INPUT INSTRUCTION 


Potential Flow Program Input 


Card 1 —Title card (required) 


Format (18A4) 
cc 1-68 TITLE Title of the case, 68 alphanumeric characters. 


cc 69-72 TITLE(18) If two points (2 sets of x,y,z) per card input, 
leave blank. 


If 1 point (| set of x,y,z) card input, punch 
"bONE", where b means blank. 


Card 2 —Control parameter card (required) 

Format (A4, 1313, 3F4.0, F12.0) 

cc 1-4 CASE Case identification 

cc 5-7 LIFSEC Total number of lifting sections 
cc MOMENT Moment origin input flag 

cc KUTTA Kutta point input flag 

cc NOFF Off-body point input flag 

cc LIST Partial execution flag 

cc MPR Matrix print flag 

GG IQUT Geometry quantities print flag 
cc IG Ignored element flag 

cc LASWAK Special last wake calculation flag 
cc IATACK Number of angles of attack 

cc IWIDTH Piecewise vorticity option flag 
cc ISAVE Stepwise pressure Kutta option 
cc IPCV Parabolic vorticity option flag 
cc SYM1 Flag for one plane symmetry 

cc SYM2 Flag for two planes of symmetry 


cc BL Flag for boundary-layer calculation 


4 


cc 56-58 


Note: 


Card 3 — Angle of attack card (required) 


Format 


cc 1-10 


cc 11-20 


cc 21-30 


cc 31-40 
cc 41-50 
cc 51-60 


Note: 


IN Unit number that contains the input data. 


For the calcuiation that this part of the report describes, 
the variables of this card should have the following values: 


LIFSEC 


] 
IATACK 


] 


BL = 1 for initial potential-flow calculation. 


BL 


0 for final potential-flow calculation 
after the new coordinates are generated. 


(6E10.0) 

ALPHAX (1) The x-component of the first input angle of 
attack. 

ALPHAY (1) The y-component of the first input angle of | 
attack. 

ALPHAZ(1) The z-component of the first input angle of 
attack. 

ALPHAX The x,y,Z components of the second input angle 

ALPHAY of attack, if any. If there is no second input 

ALPHAZ angle of attack, leave these fields blank. 


Repeat this card until total input angle of attack = IATACK 
with two angles of attack on each card. 


Card 4 — Moment origin card (optional) 


Format 

cc 1-10 
cc 11-20 
cc 21-30 


Note: 


Format 


cc 1-4 


cc 5-8 


(6E10.0) 

QMRIGNX 

PRIGNY The x,y,Z components of the input moment origin. 
QMRIGNZ 


This card is used if M@MEN40 on card 2. 


Card 5 —Lifting section information card (required) 


(1614) 


NS@RCE(1) Number of on-body elements on each strip in the 
first lifting section. 


NWAKE(1) Number of wake elementson each strip in the first 
lifting section. 
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cc 9-12 
cc 13-16 


Note: 


. aia | 


NSTRIP(1) Number of strips in the first lifting section. 
IXFLAG(1) Extra strip flag in the first lifting section. 


The same format is repeated from cc 17-cc 64. Four lifting 
sections information can be put on 1 card and may repeat the 
same format on the next card if there are more than four 
sections. Users are advised to use only 1 lifting section if 
streamline calculation is used. 


Card 6 — Ignored Element Card (optional) 


Format 


cc 1-4 


cc 5-8 


Note: 


Card 7 — Strip widths card (required) 


Format 


cc 1-10 


cc 11-20 


cc 21-30 


cc 31-40 


Note: 


Card 8 — Body definition card (required) 


Format 


Format 


3 

(1614) 

1G1(1,1) First element number on the first strip of the 
first lifting section to be ignored. 

IGN(1,1) Last element number on the first strip of the i 

first lifting section to be ignored. 


This card is needed only if IG#0 on card 2. Repeat the 
same format for the other strips of the same section. Start 
a new card for another section. 


; 
4 
: 
' 


(6E10.0) 


WIDXTR Width of the first extra strip of the first 
section, if any. If no extra strip input, 
leave blank. 


WIDTH(1,1) Width of the first strip of the first lifting 
section. 


WIDTH(2,1) Width of the second strip of the first lifting 
section. 


WIDTH(3,1) Width of the third strip of the first lifting 
section. 


Repeat the same format with 6 strip widths per card. After 
the last strip width, it should follow by the strip width of 
the last extra strip of the section whether there is extra 
strip input or not. 0.0 is used if there is no input last 
extra strip. 


(2(3E10.0,211) — For 2 input points per card 
(3E10.0, 211) —For 1 input point per card 


Y X,¥,Z components of the input point. 
Z 


STATUS Status flag 

Q = same n-line 
new n-line 
new section 
end input 


Element type 
0 = nonlifting 
1 = lifting 


Note: Repeat the same group from cc 33 - cc 64 for the next input 
point, if two points per card format were used. 


Card 9 — Kutta control point card (optional) 


Format (6£10.0) 


cc 1-10 CUTTAX( 
cc 11-20 CUTTAY ( 
cc 21-30 CUTTAZ( 
( 
( 
( 


I) 
1) X,¥Z components of the first input Kutta 
I) control points. 


cc 31-40 CUTTAX 
cc 41-50 CUTTAY 


rt} X,¥sZ components of the next input Kutta 
cc 51-60 CUTTAZ(I+1 


) control points. 


Note: This card is used only when KUTTA#0 on card 2. 
Card 10 — Kutta normal card (optional) 
Format (6E10.0) 


cc 1-10 — CNX(T) X,y,Z components of the input Kutta control 


11-20 CNY(T) point normal vector, 


21-30 CNZ(1) 


X,¥,Z components of the next input Kutta control 


pe CNY (141 point normal vector. 


31-40 any} 

51-60 CNZ(I+1) 
This card is needed only if KUTTA#0 on card 2. 

Card 11 — Off-body point card (optional) 

Format (2(3E10.0,11) 

cc 1-10 X,¥sZ Components of the input off-body points 

ce 11-20 »¥Z comp p y points. 

cc 21-30 


cc 3] LABEL Status flag 
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cc 32-4] X@FF(I+1) 
cc 42-51 YOFF (1+1) X,¥»Z Components of the next off-body point. 
cc 52-61 ZPFF(I+1) 


cc 62 LABL Status flag 
Note: LABL=3 at the end of the off-body point input, otherwise 
leave it blank. This card is needed only when NOFFA0. 
ms “se Program Input 
- .itle and option card (required) 

Foruiat (18A4,12) 

cc 1-72 TITLE Any appropriate title for the case to be run. 
cc 73-74 IMPT Option flag 

IPPT=0 No option used. 

IMPT=1 No streamline calculation but uses the 
previous streamline results to prepare 
for the small crossflow boundary-layer 
calculation. 

IPPT=2 Streamline calculation only, 

Note: When the flow-field data of the case is loaded first time 


onto the direct-access unit, set TITLE(18)="bbb1" this 
enables the program to print the data for user to check. 


Card 2 — Master directory title card (required) 

Format (111,19X,10A4) 

cer || MFLAG MFLAG=0 to initialize the master directory, 
i.e., anytime user wishes to override 
the direct-access data or defines a 
new data file. 

MFLAG4O Implies the data set has been previously 

stored. 

cc 21-60 TITLEM Any description for the master directory. 


Card 3 — Set directory title card (required) 


Format (211 ,1X,11,F6.0,10X,10A4) 
cc ot LASTS This flag should be set = 1. 
cc 2 NEWS NEWS=0 New data set input. 
NEWS#0 Data set already stored. 
ce 4 NSET Data set number. This should be set = 1 
cc 5-10 MACH Mach number 
cc 21-60 TITLES Any description of the set title. 


16 


Card 4 — Region set control card (required) 


Format (211,12 ,2F6.0,4X,10A4) 
cc 1 LASTAB LASTAB=0 Indicates more than one region. 
LASTAB#0 Means this is the last region, 
ce 2 NEWAB NEWAB=0 For new region. 
NEWAB#40 for subsequent regions. 
cc 3-4 IAB Region number 
a foie oe 3 a These values should be set =0 when input. 
cc 21-60 TITLEA Region title 
Note: Lower surface of the wina is considered as Region #1, and 


upper surface is considered as Region #2 in this program. 
Card 5 — Subregion set control card (required) 
Format (211,12,311,13X,10A4) 
cc. | LASTR LASTR=0 Indicates more than 1 subregion within 
the region. 
LASTR#0 Implies this is the last subregion. 
6G 2 NEWR NEWR=0 For new subregion. 


NEWR#0 For subsequent subregions within the 
region data set. 


cc 3-4 IREG Subregion set number. 

ce 5 IDTYP(1) This field is always set equal to 2. 

ce 6 ISQRCE This field is not used at present and is set 
to 0. 

ce 7 IRW Data read or write flag 


IRW=0 Data will be input and saved on direct- 
access unit 10. 

IRW=1 Pata has been saved and will be read in 
from direct-access unit 10. 


cc 21-60 TITLER Subregion title 


Note: Cards numbered 4, 5 used as a group for the lower surface 
of a wing. The same group of cards is used for the upper 
surface before the next card is input. 
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Card 6 — Streamline input card #1 (required) 


Format 


cc 1 


cc 2-3 


cc 4-5 


cc 6-7 


cc 11-30 


ce 31-35 


cc 36-45 


cc 46-55 


cc 56-65 


(11,312,311 ,1012 611 ,3F10.7) 


LASTR 


NDSET 


IABSET 
IR 


IN@RM 
ISURF 
IPF 
ISR(1) 
to 
ISR(10) 
ISF(1) 
to 
ISF(5) 


DELTA1 


EPS] 


EPS2 


| 
| 


Last flow region flag 

LASTR=0 This is not the last streamline set of 
data. After all the streamlines for 
this region are calculated, subroutine 
STREAM will read another region data 
from unit 10 and input. 

LASTR=1 This is the last region. After the 
streamline for this region is completed, 
no further region data will be read. 


Data set number where surface flow data will be 

found on #10. 

Note: This flag should be the same as NSET on 
card 3. 


Region set number. This number should be agree- 
able to IAB in card 4. 


Subregion number. This number should be the 
same as IREG on card 5. 


This field is set equal to 0. 
This field is set equal to 0. 
This field is set equal to 0. 
Subregion flow number. Since this version of 


the program does not use subregion, ISR(1)=l, 
and the rest of the array is set equal to 0. 


Secondary flow number. The whole array is set 
equal to 0. 


A criteria used with EPS1 to ensure that the 
streamline point has reached the stagnation line. 
Suggested value = 0.999. 


Suggested value for this criteria is 0.01. 


A criteria for terminating the streamline calcu- 
lation, if the difference between two success- 
ive coordinates is less than EPS2. Sugaested 
value 0.001. 


Card 7 — Streamline input card #2 (required) 


Format 


ce 1-2 


cc 3-4 


cc 21-60 


NSTR 


IRSAVE 


TITLER 


(212,16X ,10A4) 


Streamline save flag 

IRSAVE=0 Do not save. 

IRSAVEZ0 Streamline data will be saved on 
unit 10. 


Number of streamlines to be calculated for the 
input region number. 


Title for streamline flow region. 


Card 8 — Streamline input card #3 (required) 


Format 


cc 1-2 


cc 3-4 


ce 5 


cc 31-40 
cc 41-50 


(212,311 ,12,14,7X,4F10.0) 


IPRINT 


ISAVE 


ISTART 


ISTAG 
ISM@DE 
IPANL 


L 


DELTAS 


| 


XSI 
YSI 


Streamline print flag. | 

IPRINT=0 No streamline data will be printed. 

IPRINT=2 Print streamline data for every second 
DELTAS point. 

IPRINT=3 Print streamline data for every third 
DELTAS point. 

etc. 


Streamline save flag. 

ISAVE=0 No streamline data will be used. 

ISAVE=2 Save streamline data for every 
second DELTAS point. 

ISAVE=3 Save streamline data for every 
third point. 

etc. 


Streamline starting condition flag. 

ISTART=1 Starting point of streamline. Caicu- 
lation is obtained from unit #10. 

ISTART41 User will input the starting point 
of the streamline calculation. 


These fields should be left blank at present. 


Streamline integration distance step interval 
for the Runge-Kutta integration process. (May 
be negative to integrate forward. ) 


X,;¥,Z components of the streamline starting 
point. Input these fields only if ISTART#1. 


Repeat card no. 9 as many times as the number of streamline NSTR input in 
card no. 8. 


3.3 


Note: 


Cards nos. 7, 8, 9 are input as a group for one region (say, 
lower surface of the wing). it LASTR=0 is input in card no. 
7, then the same cards (i.e., 7, 8, 9) should be used again 
for additional region. 


Card 9 — Reference velocity and Reynolds number card (required) 


Format 


(3F15.6) 
VREF 


REREF 
SF 


Reference velocity in ft/sec to be used in 
boundary-layer program. 


Reynolds number in unit of eto) 


Input scale factor to convert the basic points 

in the unit of ft. For example, if the original 
input to the potential-flow program is in inches, 
then the SF here should be input as 12 so the 
program can convert the points to feet. SF is 
input as 1.0 if original input is already in 
feet. 


Small Crossflow Boundary-Layer Program Input 


Card 1 — Option card (required) 


Format 


ec 1-5 


cc 6-10 


cc 11-15 
cc 16-20 


ce 21-25 


(615) 
NCASE 


INTAPE 


IGUTAP 
IPT 


JQPT 


Total number of cases input. 
Example: If there are six strcamiises in one 
region (say, lower surface), then NCASE=6, etc. 


Unit number that contains the input da .u. 
Unit number that will contain the outrut. 


Print option flag. 

IMPT=0 Short print 

IPPT=1 Selected output print 
IQPT=2 Longa print 


Selects input region number and streamline 

number to be executed. This option is designed 

for users who wish to execute only the selected 

streamlines in the boundary-layer program and 

skip the rest of the streamlines. 

J@PT=0 Execute every streamline. 

J@PT#0 Input the region number and streamline 
number. 

Notice that if only 3 streamlines were selected, 

then NCASE=3 is in cc 5. 
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cc 26-30 N@PT N@PT=0 Program will use the input transition 
location. 

N@PT#0 Program will check the input transition 
location. If the input transition loca- 
tion is beyond the acceptable range, the 
program will override it with a location 
that coincides with the change of veloc- 
ity profile. 


Card 2 — Region and streamline number input card (optional) 


Format (20A4) 

cc 1-4 SET(1) Region number in which the streamline will be 
used, 

cc 5-8 SLINE(1) Streamline number to be executed in the 


boundary-layer program. 
Repeat the same format for the next selected streamline to be executed. 
The streamline numbers may be entered on one card. This card is not 
needed if J@PT=0 is entered on card 1. 
Card 3 — Region title card (optional) 
Format (16A4) 
cc 1-40 TITLE(1) 
to Any suitable header. 
TITLE(10) 
cc 41-44 TITLE(11) This field should contain region number. 
cc 60-64 TITLE(16) This field should contain streamline number. 


Note: If the data come from a prestored unit, i.e. INTAPE#5, this 
card is not required. 


Card 4 — Control flag card (optional) 


Format (511,313,3E14.6) 

cc 1 FL1 This flag is set equal to 1. 
ce. 2 FEZ This flag is set equal to 2. 
cc 63 FES 

cc 4 FL4 These flags are set equal to 0. 
cé 5 FL5 

cc 6-8 NXT Total number of points. 

cc 9-11 NTR Transition location. 


2] 


est eee ORE Os 


cc 12-14 NNN (Not used) 

cc 15-28 VGP VGP=1.14 

ce 29-42 DETA] DETA1=0.01 

cc 43-56 EPS EPS=0.001 

Note: This card is needed only if INTAPE=5. 


Card 5 — Velocity and Reynolds number card (optional) 


Format (4E14.6) 

cc 1-14 VREF Reference velocity in ft/sec. 
cc 15-28 — REREF Reynolds number in ft"! 

cc 29-42 SWANG This field is input with 0.0 
cc 43-56 P30 P30=1.0 for user input . 
Note: This card is required only if INTAPE=5. 


Card 6 — Station data card (optional) 


Format (414.6) 

cc 1-14 X(I) Streamline distance between point I and the 
starting point. 

cc 15-28 UE(I) Surface velocity at Ith station. 

cc 29-42 RK1 (i) Orthogonal curvature at Ith station. 

cc 43-56 RK2(1) Streamline curvature at Ith station. 

Note: One station data per card. Repeat the same format for 


another station. This card is needed only if INTAPE=5. 
Card 7 —Transition location card (required) 
Format (15) 
cc 1-5 NTR Streamline point number designated as the user 


input transition location. Note that point 
No. 1 is the stagnation point. 


4 Back Interpolation Program Input 


Card 1 — Data set title card (required) 
Format (16A4,14) 


cc 1-64 PETCE Header of the region to be executed. 
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cc 65-68 ISET ISET=1 For one region. 
ISET=2 For two regions. 


Card 2 — External unit definition card (required) 
Format (1415) 


ce 1-5 IQUTAP Unit number from boundary-layer program which 
Stores the 6* values. 


Scratch units to be used in the matrix 
solution subroutines. 


Y 


ce 11-15 IUNITC 


ce 6-10 IUNITB 
cc 16-20 IUNITD 


cc 21-25 IFT20 Unit number from the potential-flow program which 
stores the Ug, Vo values. This unit number is 
defined as no. 20. User may change it in the 
potential-flow program if so desired. 


cc 26-30 IFT22 "nit number that stores the interpolated 6* at 
the potential-flow calculation control points. 
This is the output unit from this program and it 
will be used in the final potential-flow 
computation. 


cc 31-35 KS@RCE Number of source elements per lifting strip as 
input to the potential-fiow program. 


cc 36-40 KWAKE Number of wake elements per lifting strip as 
input to the potential-flow program. 


ce 41-45 ITEXTRA Number of extra strips as input to the potential 
flow program. 


Card 3 — Streamline information card (optional) 
Format (1415) 
ce 4-5 ISTRM(I), Streamline region number, 


ce 6-10 NSLINE(I) Number of streamlines within the region, 


cc 11-15 LSTART (I) Starting address of the first streamline within 
the region, 


cc 16-20 ILIM1(1) Starting element number of the strip from | 
potential-flow program that is used as the lower | 
limit of the region. 


ce 21-25 ILIM2(1) Ending element number of the strip from the 


potential-flow program that is used as the upper 
limit of the region. 
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hatte ee ee 


This card is needed only if user wants to run this program as a 
separate run. Otherwise data will be stored and passed into 
this program. 


Card 4 —u,v data card (optional) 
Format (6F12.6) 


cc 1-12 UTE(I) u-coordinate of the u-v plane at the trailing 
edge. This is calculated in the potential-flow 
program. 


ULE (I) u-coordinate at the leading edge. Obtained from 
potential-flow program. 


VTE(I) v-coordinate of the u-v plane at the trailing 
edge. This is calculated in the potential-flow 
program. 


VLE(I) v-coordinate at the leading edge. 


CROPT(1) Distance of the root chord in the u-v plane 
between the trailing edge and the leading edge. 
This is obtained from streamline proaram. 


cc 61-72 CTIP(I) Distance of the tip chord in the u-v plane 
between the trailing edge and the leading edge. 
This also is obtained from streamline program. 


Note: All the above data is for region I. This card is needed 
only if the program is run separately. 


Card 5 — Streamline output condition card (optional) 


Format (1415) 


ce 1=5 NSPT(I,IS) Number of streamline points of the Ith streamline 
in the ISth region. 


IC@ND(I,IS) ICQ@ND=0 Means the streamline computation ends 
normally and the points can be used. 
ICQ@ND=1 Indicates the streamline computation 
does not end normally and this stream- 
line should be skipped. 
I For streamline count. 
IS For region number. 


LP@INT(I,IS) Starting address in the direct-access unit #10 
to be used for retrieving the streamline data. 


One streamline information per card. Repeat the same format 
for another streamline. This card is needed only if the 
program is run independently. 


3.5 Final Potential-Flow Program Computation 


This part of input is the same as the initial potential-flow program (see 


3.1) with the following exceptions: 
1. BL on card 2 is input with a negative value, say -1.0. 
ae Body points input is not needed. 


The remaining input is the same as the initial potential-flow calculation 


input. 


4.0 DECK SETUP 


Two types of JCL deck setups are given for a multistep job. The first 
type shows the JCL setup for computing, linkediting and executing each step 
of the program. The second JCL setup is for the case where all the steps of 
the program have been stored on a diskpack prior to execution. 


Although these examples are for u multistep job, each program may be run 
individually. In this case, the data sets to be saved in each step should be 


Dati line 


kept rather than passed into the next program step. For example, in the 
potential-flow program, initial calculations, there is a DD-card for Unit 21 
defined as follows: 

//GO.FT21F001 DD DSN=RUNCASE.FL@WDATA,UNIT=TAPE16, 


ff DISP=(NEW,PASS), 
iy) DCP=(RECFM=VBS ,BLKSIZE=6447 , BUFNQ=1), 
// LABEL=(1,SL,RETPD=180) 


The disposition parameter indicates that this data set is to be passed into 
the next step, i.e., the streamline calculation program. However, if it is 
desired to save this data set and use it at a later time, this can be 
accomplished by changing the disposition parameter to: 


DISP=(NEW,KEEP), 


in the aforementioned DD-card. 


It is also possible to allow a certain data set to be punched on cards 
rather than saved on tape or disk. In that case, the punch cards may be input 
to the next program after the output of the previous program has been examined. 
This capability is applicable to those data sets that have the record format 
RECFM=FB defined in the DD statement. 


In short, if these program steps are to be run independently, it is neces- 
sary to change some JCL parameters. Users may have to alter the JCL decks 
according to the system they are using. 
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(b) The Streamline Calculation Program 


Card Columns 


0 l 2 3 : 4 5 5 7 
17345678901 234567890) 2345678991 23456 78901 23245672901 22445672901 23456789012 


JISTEP2 EXFC FORTXCL Ge PARM FOR T=#NPT=H2,XREF TP EGTON.GO=2 50K 
JIFEQRT.SYSUT2 OD DSN=ESVSUT2,UNTT=SYSDA «SPACE RHCTRK 4 (5559) 
J/FORT.SYSIN 9D * 

( @FPLACE THIS CARD WITH THF STREAMLINE CALCULATICN PROGRAM NECKS.) 
/* 
JI/GO.FTOLFOOL OD DSN=ECUTOATA.LTL 4EY,UNITH=TAPE 16 DDISO=H(NEW PASS)» 
// NCB=(RECFM=FR, LRECL = ah, RLKSTZ7E=4900)- 
4/ LABEL =(1)SL,R2£TPN2190) 
J/GO.6FTO2ZFOOL CO GSN=EFTO?,UNT THSYSDA,SPACFHITRK (50,1 ODD, 
4/ UCB=(FECFM=VRS,8LKSI7E =6 6447» BUFNQ= 1} 
J/GO*FTOZFOOL DNDN OUSN=EFETIZA,ZUNETHSYSNIA, SPACE =(TRK (50 10)), 
// DC R=( PECEM=V359 RC KSI7E 26447 ¢RUFNOZL) 
J/COCFTOSFOOL DD DSN=RETO4SUNETHASYSIA, SOACE=ITRK (50610)), 
4/ OC B=CPECEM=VRS pBLKSIZ&=6447, BUFNO=L) 
4/GOeFTOTFOOL OD SYSNUT=8 
J/GO.FTIOFOOL DD CSN=LRACLEDH3.9PMSILHB ASTINCASE »UNTT=3230-1 4 
// DCR=(OSNRG=NA, RUFNO=L) »VILESER= CSLB304 D1 SP=(NE Wy PASS) 
// SPACE=(109,(2900)) 
JIGOFTISFOOL OD DSN=*.STE°LGOLFTO2ZF O01 »DOISP=(CIN,PASS) 
J/GOsFETL6FOOL PD USN=*.STEPL.GO AFTOGFOOL «DL SP=(QLD,PASS) 
JIGOSFT2ZLFOOL ON ASNE*.STEPLLGOLFTALEOOL DI SO=(NLN KEE) 
7/G6C.SYSIN DP * 

( REPLACE THIS CARD WITH FHE INPUT DATA TO THE STREAMLINE PROG. } 
/* 
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(c) The Small Crossflow Boundary-Layer Program 
Card Columns 


fs 
123456739912345679901234567P90173% 567890 123456729012345 6789912345789012 


FURTYCLGy2ARM FORT =Y¥REF pOPTH2*,REG 
DIY NSN=ESYSUT2,UNTT=SYSNA, SPACE=( 

no * 

THIS CARD WITH THE BOUNDARY LAYER 


CO OSN=*.STE22?.50. rir oor 


ON  OSN=NSTAP EL GE PUNT T =T re 


00 
No 
Do S¥SQU 
* 


-E THIS CAPD WITH THE INPUT DATA TO THE BOUNDARY LAYER PROG) 


A 
KS 
bo | ‘ : 1 
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+ ~BKHKNNBNNNN FH WN 


/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
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(d) The Back Interpolation Program 


FOR TXCLGeP ARM CFOR T=" XREF gOPT=2* REGION. GO=220K 
CO DNSN=ESYSUT2,UNIT=SYSDA SPACE H=(TRK,(555)) 
xc 


HTS CARD WITY THE BACK INTERPOLATION PROGRAM DECKS. ) 


Dp * .STE%3.60 2FETO2ZFOOL s+ DI SP=( OLD, KEEP) 
ETO2,UNIT=SYSDA, SPACE=(TRK y(50yL0) Dy 
VES 9 BU ite 6447, RUFNO= L) 

YSDA, SPACE HCTRK (50919) Do 
7£=6447,RUFNO=L) 
Y 
Z 


nny 
pe Fay 
=~ 12ET 
A440 


N° 
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e ee 
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—a+4 + + -+4-4 Dee $ 
=-2 oO GCG OO MmMNN 


SOAs SPACE=CTRK 9 (50510) 09 
F=6447,BUFNO=1) 
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(e) The Final Calculation Using the Potential Flow Program. 


Note that the potential-flow program is assumed to be stored on a 
volume identified as "CSLB30" with the data set name 
"DAC .EOH3. DPM.J1HB.V3". 


0) l 2 3 4 5 6 7 
12345578901 234567990123456789012 3456789012 3456789012345673991 23456789012 


J/STEP5 EXEC FORTCO,REGIIN=380 
J/STEPLIR OCD OSN=DAC.EOH2.DPM.LILHB.V2gDI SO =H(SHR,PASS) FUNIT=33320—-1 + 
4/ VOL=SER=CSLBA30 
4/GC.SYSIN ODO * 
4/GC.FTOILFOOL ND OSN=KETIOLZUNTT=HSYVSDA, SPACER TRK 5 (50910)), 
4/ DCB=(RECEM=VBS, A8LKSITZE=6447,BUFNO=1 ) 
J/GO.FTO2FU0L OD OSN=EFTO2Z,UNET=SYSDA, SPACE=(TRK (50910) ),% 
Vf NCR=*.FTOLFIO) 
J/GCeFTO3FOOL DD NSN=EFTO3,UNITHSYSDA, SPACEH(TRK (50919) 04 
4/ NCRB=*. FTOLFOO! 
//GO.FTO4FOOL OD OSN=EFTO4, UNI T=SYSDA, SPACE=H(TRK 6150410) ) > 
4/ NC B=*. FTOLFIO} 
//GC.eFTOTFOOL OD SYSOUT=8 
4/GC.FTO8BFIOL ND DOSN=EFTOBSIINE T=SVYSPA,SPACFEH(TRK (50,10) D4 
// NCB=*.FTOLFOOL 
J/GO.FTO9SFODO1 NN ONSN=EETISGLUNE T=SYSNA, SPACE H=H(TRK 6 (50010) De 
// NOCB=*.FTOLFIOL 
J/GOeFTIOFIOL NY OSN=SFETLOSUNTT=SYSOA, SPACE=H(TRK 9 (50910) Dy 
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( REPLACE THIS CAPD WITH THE IN®UT DATA TO’ THE ETNAL PLFLOW PROG. ) 
/* 
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4.2 JCL Setup for Executing the Programs Already Stored on an External Unit. 


All program steps are stored in the volume "CSLB30." 


(a) The Potential-Flow Program Step has the data set name 


"DAC. EQH3.DPM.J1HB.V3." 
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HIS CARD WITH THE INPUT DATA TO POTENTIAL FLOW PROG. 
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The streamline Calculation Program has the data set name 
2 


“DAC. EOH3.DPM.J1HB.STREAM." 
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(e) The Final Calculation Using the Potential-Flow Program has the 
data set name "DAC.EOH3. OPM.J1HB.V3." 
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HIS CARD WITH THE INPUT DATA TO THE FINAL P.FLOW PROG. 
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5.0 SAMPLE TEST CASE 


The sample test case is an eight-strip swept wing with 30 on-body elements 
and one wake element per strip. The test case is run at an angle of attack of 
8.22° with one plane of symmetry. 


The input data and partial output are included. Due to the size of the 
output from the programs, only the output for the first strip is presented 
here. 
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Input Data to the Potential Flow Program. 
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(c) Input Data to the Boundary-Layer Calculation Program. 
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(d) Input Data to the Back Interpolation Program. 
NACA 8 STRIP SWEPT WING. BACK INTERPOLATION. BOTH SURFACES. 2 

L 2 3 4 20 22 30 l 0 
1 8 504 1 15 

0.0 46.67 0.0 0.0 42.1256 21.0475 
fs 8 1167 16 20 
65 0 517 
62 0 582 
63 0 644% 
Tl O 707 
69 Q 778 
66 0 847 
59 0 913 
18 1 972 

18 1 1180 
43 0 1198 
40 0 1241 
40 O 1281 
41 0 1321 
42 QO 1362 
42 0 1404 
41 QO 1446 
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Input Data to the Final Calculation in the Potertial Flow Program. 


NACA SWEPT WING, 8 STRIPS,» A=8.22¢ 30 SOURCE, 1 WAKE, L PLeSYMe PeKUTTA 
NACA | i \ { 1.0 0.0 1.0 
0.989707 O. 0.143106 
30 1 8 0 
° 2-17 2-80 4.70 6.10 10.00 
8.40 6.80 «10 0.00 


——— ee 


5.2 Output Data 


The output data from the test case consists of the output from these 
programs : 


1. The initial potential-flow calculation (pp 41 - 49). 

2. The streamline calculation (pp 50 - 68). 

3. The small cross-flow boundary-layer calculation (pp 69 - 73). 

4 The final potential-flow computation (pp 74 - 82). 

Output of the first strip from this eight-strip test case is used here 
as a sample. Similarly, the first streamline output of the lower surface of 


the wing is selected in both streamline calculation and the boundary-layer 
calculation. 
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